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Abstract 

The BABAR detector at SLAC's PEP-II storage ring has collected data equivalent to about 30.4 fb _1 
through June 2001. Results on CP violation, and in particular searches for direct CP violation, 
and measurement of rare B decays are presented. 
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1 Introduction 



We present a sample of recent results from BABAR including the observation of CP violation with 
the measurement of sin2/3 and of rare decay modes. 

The BABAR experiment has been running at the PEP-II asymmetric e + e~ collider since 1999. 
Its main goal is a high statistics study of B decays, including a study of CP violation in the B sector. 
The center-of-mass energy is tuned to the T(4S) just above the BB threshold. The asymmetry of 
the energies of the two beams provides a longitudinal boost so that the average B flight length is ~ 
250 fim and can be measured. The high nominal luminosity of 3 x 10 33 cm - 2 s allows the study 
of many rare B decay channels. 

The BABAR detector is described elsewhere [[jj. Charged particle track parameters are obtained 
from measurements in a 5-layer double-sided silicon vertex tracker and a 40-layer drift chamber 
located in a 1.5-T magnetic field; both devices provide dE/dx information. Additional charged par- 
ticle identification (PID) information is obtained from a detector of internally reflected Cherenkov 
light (DIRC) consisting of quartz bars that carry the light to a volume filled with water, and 
equipped with 10752 photomultiplier tubes. Electromagnetic showers are measured in a calorime- 
ter (EMC) consisting of 6580 CsI(Tl) crystals. An instrumented flux return (IFR), containing 
multiple layers of resistive plate chambers, provides fi identification. 

Most results presented here are obtained with a data sample of 30.4 fb _1 collected through 
June 2001, for a total amount of about 32 x 10 6 BB pairs. 

2 Measurement of sin2/3 

The primary goal of the experiment is the observation of CP violation in the B sector. In the 
standard model (SM), CP violation occurs via a complex term in the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. The phase of such a complex term can eventually be measured from the interference 
of two amplitudes contributing to the same final state. In the Wolfenstein parametrization[Q] of 
Vckm, all the matrix elements are real, except V t d and V u b. Most of the envisaged strategies to 
observe CP violation therefore use the interference between an amplitude containing one of these 
target complex CKM elements and a real amplitude, to a CP eigenstate. 

In the case of a b — > ccs decay like B -> J/^K Q S , the B can either decay directly via a real 
amplitude oc V* b V cs , or decay after a B° — > B° oscillation that proceeds through a box diagram 
with an amplitude (V^jV^) 2 with a phase —2/3. The measurement is theoretically clean, as the 
main higher-order diagrams have the same weak phase as the first-order ones. It is also relatively 
background free due to the presence of a J/ip in the final state, and benefits from a branching ratio 
,8~ 10~ 3 , a rather high value in the land of B decays. 

The resulting evolution function g±(t), for a decay to a final state /, involves a complex quantity 
Xf. If the (so called direct) CP violation in the decay itself can be neglected - as it is the case here, 
within 1 %, in the standard model - we have \Xj\ = 1, and we obtain: 

g±{t) = x [1 ± ImXf sm(Anidt)] , 

TB 

where t is the B proper decay time, and the + or — corresponds to an intial B° or B° meson. 
In the present case of a b — > ccs decay, Xf is 7]fe~ 2t ^, where r\f is the CP eigenvalue of /, and g 
simplifies to: 

g±(t) = x [1 T ??/sin2/3sin(Am d i)] , (1) 



2 



allowing the measurement of sin2/3. The "initial" state of the CP B can be known thanks to the 
EPR paradox: the T(AS) resonance having spin-parity 1 , decays to a coherent BB pair in an 
antisymmetric state, so that when the first B decays, say to a B°, the other one is in the opposite 
state, a B°. The oscillation time is then the difference At of the time of flight of the two B's, 
measured from the difference Az of their flight lengths. The vertexing of the CP B is performed 
with an excellent precision, a z ~ 60 /im, thanks to the presence of the hard leptons from the 
J/tp . The other B is vertexed inclusively, with a x 2 cut that limits the systematics due to cascade 
charmed decays; the contribution of the other B dominates the resolution on Az, which is of the 
order of 180 /um independent of the CP channel used. 

The flavor determination of the other B, called tagging, is also performed inclusively: events are 
classified in four mutually-exclusive categories, using respectively the charge of the fastest lepton, 
the total charge of identified kaons, the output of neural networks that use the information carried 
by non-identified leptons and kaons, and by soft pions from D*'s, or are not tagged ||. The effective 
tagging efficiency Q{ of each category i is defined as the product of its efficiency ei and of the square 
of its dilution (1 — 2i«j). The total effective tagging efficiency Q = J2i=i e «(l ~~ 2u>j) 2 is measured on 
a large sample of exclusively reconstructed B decays of specific flavor (B — > D^~h + , h = ir, p, a\, 
and B -> J/0 K*°(K + tt-)), and is equal to (26.1 ± 1.2) %. 

With a real detector, sin2/3 is multiplied by a dilution factor (1 — 2w) in eq. |l], and g±(At) 
is convoluted with a At resolution function. The value of sin2/3 is obtained from an unbinned 
maximum likelihood fit to the At distribution of a combined event sample consisting of the CP 
sample and of the flavor sample. (The number of tagged events, the purity, and the CP effective 
value of the channels used are given in table |). In this way, the mistagging probabilites ui, and 
the parameters of the resolution function are determined from the data, and their correlation with 
sin2/3 is taken into account. 

Background events are taken into account by a separate pdf that enters the likelihood. When 
the CP B has no K® in the final state, the probability for an event to be a signal event is computed 
from the value of the energy substitued mass tues = ( s /4 — Pb) 1 ^ 2 where s is the total energy 
and p* B the measured momentum of the B candidate, in the T(AS) rest frame (fig. |l|a). For a CP 
B with a K® in the final state, only the direction of the K® is measured in the detector, and a 
kinematic fit of tties to the nominal B mass is performed. The signal probability of the event is 
computed from the difference AE = E* B — E£ earn between the energy of the B candidate and the 
beam energy, in the T(4S) rest frame (fig. |l|b). 

Table 1: Number of tagged events, signal purity and CP content of the CP and flavor samples. 





Channel 


Ntag 


Purity (%) 


(Vf) 


CP sample 


J/0 K% (vr+ 7T - ) 


316 


98. 


-1 




J/0 K° s (tt° vr°) 


64 


94. 


-1 




0(25) K° s (tt+ tt-) 


67 


98. 


-1 




Xcl K° s (7T+ TT") 


33 


97. 


-1 




J/0 K*° (ifO 7T°) 


50 


74. 


0.65 ± 0.07 




J/0 K° L 


273 


51. 


+1 


non CP sample 


Flavor 


7591 


86 





The slight asymmetry of the time distributions (fig. |2|) is barely visible, but the effect is clear 
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Figure 1: a) Distribution of tties f° r B candidates having a K® in the final state; b) distribution 
of AE for J/ipK® candidates. 



in the plot of the asymmetry itself Acp(At) = —(1 — 2w)rjf sin2/3sin(A77i^At). We obtain sin2/3 = 
0.59 ± 0.14(siai) ± 0.05(syst). If there were no CP violation, i.e. (3 = 0, the probability of such a 
measurement would be 3 x 10 -5 . Therefore, our result for sin2/3 represents a 4.1 standard deviation 
observation of CP violation j|] in the neutral B meson system. 

The impact of the present measurement on the knowledge of the CKM matrix is shown in fig. 
H| The set of measurements is clearly consistent within the framework of the standard model, while 
the present measurement significantly decreases the size of the la allowed region. 

We have also searched for direct CP violation in the decay. Releasing the constraint [] |A/| = 1, 
the evolution equation becomes: 



3±(At) 



e -|At]/r fl 

2r B (l + |A2|) 



l + |A/r\ / 1 



± f-^C 1 - |A/| 2 )cos(Am d At) + lm\ f sin(Am d At)J 



Fitting for |A/|,ZmA//|A/| with the r// = —1 sample, we obtain[||] |A/| = 0.93 ± 0.09 ± 0.03, 
showing the absence of direct CP violation at the 10% level. 

3 Charmless two body B decays 

The V u b term that appears in the amplitude of the decay b — > u carries a 7 phase: this makes the 
measurement of sin2a in B — > tt + it~ decay-mixing interference attractive, using a time distribution 



but still assuming that there is no CP violation in mixing. 
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Figure 2: Left: Time distribution of rjf = — 1 candidates with a B tag N B o and with a .8° tag 
Ngo, and the asymmetry (iVgo — N^ )/(N B o + Ng ). The solid curves represent the result of the 
combined fit to all selected CP events; the shaded regions represent the background contributions. 
Right: Corresponding information for the r]t = +1 mode (J/tpK®). 



analysis similar to that for sin2/3. However, in this case the decay has a strong CKM suppression 
(V u b oc A 3 ), so that the higher-order diagrams, so called penguins, contribute at the same level as 
the tree diagram. This leads to a complication in the measurement of sin2a, but also makes the 
process sensitive to the contribution of new heavy objects (Higgs, SUSY) in the penguin diagram, 
which would result in direct CP violation via tree-penguin interference in B — > K + ir~ decays. 

On the experimental side, the difficulty is background rejection. The decaying 5's are almost 
at rest in the T(45) rest frame, and their 2-body decay products are approximately back-to-back 
in that frame, while continuum (e + e~ — > qq) events have a two-jet topology that provides a lot 
of back-to-back track pairs. BB decay products are roughly isotropically distributed in the T(4S) 
frame and we use this behaviour to separate signal from noise. 

Preselection cuts on global event shape parameters (Fox- Wolfram moment, sphericity), on the 
angle between the sphericity axes of the B and that of the rest of event, and on a Fisher discriminant 
J- that optimizes the use of the energy flow with respect to the decay axis §. Kaon/pion separation 
is achieved using the Cherenkov angle 9 C measured in the DIRC and the value of AE (computed 
with the pion hypothesis for both tracks). 

The branching fractions are obtained by an unbinned maximum likelihood fit with the 4 
variables: m-ES) AE, J 7 , 6 C and are given in table [| The corresponding ttjes and AE spec- 
tra, after appropriate cut on the likelihood ratio to enhance the signal fraction, are given in 



5 



-1 1 2-1 1 2 

P P 

Figure 3: Allowed region (light (dark) grey: two (one) standard deviation) in the CKM parameters 
p, r] plane. Left: the present sin2/3 measurement is overlaid on preceding fit. Right: the present 
sin2/3 measurement is included in the fit. The fit method is described in ref. ||. 



Table 2: Detection efficiencies (e), fitted signal yields (Ns), statistical significances (S), branching 
fractions (B), and charge asymmetries of the charmless two body decays (20.7 fb _1 ). 



Mode 


e(%) 


N s 


S(a) 


B/W- b 


A C p 




45 


41 ± 10 ± 7 


4.7 


4.1 ± 1.0 ±0.7 




K+7T- 


45 


169 ±17 ±13 


15.8 


16.7 ±1.6 ±1.3 


-0.19 ±0.10 ±0.03 


K+K- 


43 


8.2j£! ± 3.5 


1.3 


< 2.5 (90% C.L.) 




7T+1T 


32 


37 ± 14 ± 6 


3.4 


< 9.6 (90% C.L.) 




K+TT° 


31 


75 ± 14 ± 7 


8.0 


10.8tH ± L0 


0.00 ±0.18 ±0.04 


K°TT+ 


14 


59±\o ± 6 


9.8 


I8.2±fl ± 2.0 


-0.21 ±0.18 ±0.03 


K°K+ 


14 


-4.1±f;| ± 2.3 




< 2.4 (90% C.L.) 




K°ir° 


10 


I7.9±f;l ± 1.9 


4.5 


8.2jJ;£ ± 1.2 





fig. ||. Self-tagging modes are used to search for direct CP violation in the charge asymmetry 
A CP = (B(B - /) - B(B -> f))/(B(B - /) + B(B -» /)). 

A vigorous effort is under way to understand the shift in sin2a due to the presence of the 
penguins 0. The quantity that is actually measured is denoted sin2a e ff. 

An unbinned maximum likelihood fit similar to that used for the branching fractions is used 
to measure sin2a e jj: the variable At is added in the likelihood. The vertexing and tagging are 
similar to that used for sin2/3. There is here a (strong) probability of direct CP violation so that 
the evolution equation used is: 

e -\At\/r B 

g±(t) = — x [1 ± (-C f cos(Am d t) + S f sm(Am d t))] 

Fitting for the coefficients Cf = (1 - |A/[ 2 )/(l + |A/| 2 ) and S f = 21m\ f /(l + \X f \ 2 ) we get g, 
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Figure 4: The rriES and AE" distributions of charmless two body decays, after a likelihood ratio 
cut. The solid curves represent the fit predictions for both signal and background; the dashed curve 
represents the given signal mode only and the dotted curve represents other modes of the same 
topology (20.7 ftT 1 ). 



on 30.4 fb _ \ S f = 0.03±g;|g ± 0.11 and C f = -0.25±[J;^ ± 0.14 (preliminary), where S f would be 
equal to sm2a e ff if |A/| = 1 (No direct CP violation). 

4 Radiative penguin decays 

b — > decays like B — > K*j are of particular interest because flavor changing neutral currents 
(FCNC) are forbidden in the standard model. Effective FCNC are induced by one-loop penguins 
in which the top quark dominates. Again, unknown heavy objects like charged Higgses or SUSY 
partners can contribute in the loop, and the interference with the SM diagram can induce CP 
violating charge asymmetries as large as 20 %, while the SM predicts Acp < 1 % |J. 

Monochromatic photons are selected (2.30 < E* < 2.85 GeV). As already noticed before, a 
large background from continuum events is present in this two body decay mode. The direction of 
the photon and the thrust of the rest of the event, in the T(4S) rest frame, are uncorrelated for 
B — > K*^f events, and strongly correlated for background, and a cut on their respective angle is 
the most powerful means of rejection available here. The branching fractions are measured with 
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Table 3: Number of events, significance, and branching fraction of B — > K*j decays (20.7 fb ). 



Mode 


B{B -> K*7)/10- 5 


K+ TT- 


4.39 ±0.41 ±0.27 


K° s vr° 


4.10 ± 1.71 ±0.42 


K° s 7T+ 


3.12 ±0.76 ±0.21 


K+ 7T° 


5.52 ± 1.07 ±0.33 




Figure 5: uies distributions of the four K*j channels (20.7 fb ). 



Charge asymmetry is measured for self-tagging modes (i.e. not 7X3^°) and is found compatible 
with zero, Acp = —0.035 ± 0.076 ± 0.012: this measurement already constrains possible physics 
beyond the SM. 

5 Gluonic penguins 

The B — * <pK^ decays present several interesting aspects: the decay is dominated by a gluonic 
penguin contribution, making this mode a smoking gun for the observation of this processes; the 
CP final eigenstate and the real decay amplitude make the measurement of sin2/3 possible, in a 
system that has a different sensitivity to physics beyond the SM compared to b — > ccs decays [jHJ . 
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Table 4: B — > (frK^ fitted number of signal event n^g, statistical significance S, and measured 
branching ratio £>. The subscripts in the <pK* + modes refer to the kaon daughter of the (j)K* + (20.7 
fb- 1 ). 



Mode 


n sig 


5 


£(io- d ) 


<f>K+ 
4>K° 




10.5 
6.4 


7.7t^ ± 0.8 
8.1±|;| ± 0.8 


<pK*+ 




4.5 


9.7+fJ ± 1.7 


<^ + + 
4>K*+ 


4-4±i:S 


2.7 
3.6 


12.81^ ± 3.2 
8.0+1? ± 1.3 


(j)K* 




7.5 


8.7^i ± 1.1 



Branching fractions are measured (table |4j) with an unbinned maximum likelihood fit to mES) 
A.E, m(K + K~) (and i^7r mass and K* helicity angle for <j>K* modes). The decays B + — » (/>.fr* + 
and -B — > i^i^ are first observations] 10]. 



6 Summary 

We have presented a sample of recent results of the BABAR experiment. CP violation is observed 
in the B sector with a significance of 4.1 standard deviations. The observation of several rare B 
decays is reported, several of which allowed a search for direct CP violation in charge asymmetries. 

The measurements presented here are statistically limited. The increase in integrated lumi- 
nosity that is expected over the next few years - 500 fb -1 in year 2005 - will allow a significant 
improvement in the precision of these results. 
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